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The effective action associated with the trace anomaly provides a general algorithm for 
approximating the expectation value of the stress tensor of conformal matter fields in arbi- 
trary curved spacetimes. In static, spherically symmetric spacetimes, the algorithm involves 
solving a fourth order linear differential equation in the radial coordinate r for the two scalar 
auxiliary fields appearing in the anomaly action, and its corresponding stress tensor. By 
appropriate choice of the homogeneous solutions of the auxiliary field equations, we show 
that it is possible to obtain finite stress tensors on all Reissner-Nordstrom event horizons, 
including the extreme Q — M case. We compare these finite results to previous analytic 
approximation methods, which yield invariably an infinite stress-energy on charged black 
hole horizons, as well as with detailed numerical calculations that indicate the contrary. The 
approximation scheme based on the auxiliary field effective action reproduces all physically 
allowed behaviors of the quantum stress tensor, in a variety of quantum states, for fields of 
any spin, in the vicinity of the entire family (0 < Q < M) of RN horizons. 
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I. INTRODUCTION 

The evaluation of the energy-momentum-stress tensor T*^ of quantum matter in curved space- 
times is important for understanding the possible backreaction effects of matter on the large scale 
geometry of spacetime. Quantitative control of the stress tensor is needed especially in black hole 
and cosmological spacetimes with event horizons, where general considerations indicate that vac- 
uum polarization and particle creation may lead to significant quantum effects which are cumulative 
with time. Such secular, macroscopic effects of quantum matter on the geometry of spacetime pro- 
vide the intriguing possibility of observable consequences of quantum gravity at accessible energy 
scales, far below the Planck scale. 

In the direct method of computing the quantum expectation value of the stress tensor, the 
first step is to solve the appropriate matter field equations, for a complete set of normal modes. 
With these solutions in hand, the Fock space of the quantum theory is constructed, a specific 
"vacuum" state |^') in the Fock space chosen, and the expectation value (^|ri,|\I') evaluated in 
the selected state, component by component, as a sum over the normal modes of the field. Since 
T*^ is a dimension four operator in four spacetime dimensions, the mode sum for its expectation 
value is quartically divergent. Hence a delicate regularization procedure, such as point-splitting 
must be employed in order to identify and remove the short distance divergences in the mode sum, 
absorbing them into appropriate counterterms up to dimension four in the gravitational effective 
action [1]. Only after this regularization and subtraction procedure is performed can finite results 
for the renormalized (T^j^) of physical interest be extracted. 

Since the wave equation, mode functions, and stress tensors are different for fields with different 
spin, this procedure must be carried out independently for each quantum field of interest. Likewise, 
if one wishes to consider different quantum states in the Fock space, with different boundary 
conditions on the mode functions, the calculation must be repeated for each state. Because of the 
intricacy of the subtraction procedure, together with the numerical solution of the mode equations, 
which is usually required, it is often difficult to anticipate the general physical features of the result. 
Finally, if the geometry is modified, or allowed to respond dynamically, the entire calculation of 
(T"^) would have to be repeated for each new geometry and/or at each new time step. This direct 
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method of calculating {T'\) is thus both time and computation intensive, and has limited the 
number of results for the stress tensor in fixed backgrounds to only a handful of special cases, 
making the consideration of dynamical black hole spacetimes varying with time in response to (T'^) 
prohibitively difficult, even in the case of exact spherical symmetry. 

Because of the difficulties involved in direct evaluation methods, considerable interest attaches 
to developing general and reliable algorithms for approximating (T'^) in general curved spacetimes. 
An approximate method that can capture the secular, macroscopic quantum effects on the geometry 
would be particularly interesting for applications to both cosmological and black hole spacetimes. 

In a previous article a general approximation scheme based on the effective action and stress 
tensor obtained from the conformal or trace anomaly was introduced [2]. Although the effective 
action associated with the trace anomaly is not unique, as it is defined only up to arbitrary con- 
formally invariant terms, a minimal generally covariant action can be found by direct integration 
of the anomaly [3j, The logarithmic scaling behavior of the effective action associated with 
the anomaly separates it from any of the other possible local or non-local Weyl invariant terms 
in the exact effective action, which do not share this logarithmic scaling property. The non-local 
anomaly action may be cast into a local form in a standard way by the introduction of one or 
more scalar auxiliary degree(s) of freedom 

aaaa. 

Since there are two distinct cocycles in the 
non-trivial cohomology of the Weyl group in four dimensions [3], the most general representation 
of the anomaly action is in terms of two auxiliary scalar degrees of freedom, each satisfying fourth 
order linear differential equations of motion (j2.8p . These are two new scalar degrees of freedom 
in the low energy effective theory of gravity not present in the classical Einstein theory. Since 
the effective action expressed in terms of the auxiliary fields is a spacetime scalar, variation with 
respect to the metric yields a covariantly conserved stress tensor, which not only reproduces the 
trace anomaly but also yields non-trivial tracefree components as well. 

In the auxiliary field approach, computation of the quantum expectation value {T'\) is reduced 
to the solution of linear, classical equations for the auxiliary fields, bypassing completely any 
summation over modes, and the regularization and renormalization that requires. Different states 
of the underlying quantum field(s) are associated with the choice of specific homogeneous solutions 
to the linear differential equations satisfied by the auxiliary fields. This allows for states obeying 
different boundary conditions on the horizon to be studied simultaneously. In addition, matter 
fields of every spin are treated in a unified manner, since the auxiliary field equations do not 
depend on the spin of the underlying quantum field. The stress tensor depends on the spin of the 
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fields only through the known spin dependence of the two numerical coefficients appearing in the 
anomaly in eqs. ()2.3p below. Finally, the auxiliary field action and the stress tensor derived from 
it can be evaluated in principle in any spacetime, dynamical or not, without respect to special 
symmetries. Thus the scalar effective action of the anomaly furnishes a general classical algorithm 
for approximating the expectation value of the full stress tensor of quantum matter of any spin in 
an arbitrary curved spacetime. 

The auxiliary fields are sensitive to macroscopic boundary conditions and the presence of causal 
horizons, so it is particularly interesting to apply the approximation algorithm based on them to 
spacetimes with event horizons, where quantum fluctuations are expected to play an important 
role, and comparison with existing numerical results is possible. Even when the quantum matter 
fields are not strictly massless, their fluctuations and stress tensor in the vicinity of an event 
horizon can exhibit conformal behavior [2|]. In particular, quantum states with diverging (T'^) on 
the Schwarzschild horizon, for which the backreaction on the classical geometry is significant, have 
precisely those diverging behaviors prescribed by the stress tensor derived from the effective action 
of the anomaly. 

In Ref. 0] a study of the stress tensor obtained from the anomaly in Schwarzschild spacetime 
was undertaken. In Ref. [3] the general form of the stress tensor due to the conformal anomaly in 
an arbitrary spacetime was given and applied to a few special cases, such as Schwarzschild and de 
Sitter spacetimes. A detailed comparison of the two studies is given in Section 4. 

In the present article we extend and develop the classical approximation technique of Ref. for 
the quantum stress tensor based on the trace anomaly to static, spherically symmetric spacetimes, 
focusing specifically on electrically charged Reissner-Nordstrom (RN) black hole spacetimes, and 
states with regular stress tensors on the RN horizon. Since direct computations of the renormalized 
stress tensor expectation value {TJ') have been carried out in both Schwarzschild and Reissner- 



a la, 



3, 11, 12, 3, 14, 3, i6,Q,Q, 



we 



Nordstrom (RN) spacetimes for free fields of various spin 
will be able to compare the results of the new approximation scheme to these direct computations 
of (r^^). Since the exact effective action of quantum matter generally contains terms which are not 
determined by the anomaly, using it to compute the stress tensor expectation value is certainly 
an approximation, which will differ from the direct evaluation of (T'^) in terms of mode sums 
in general. It is therefore a non-trivial check of the approximation scheme if all of the allowed 
behaviors of the traceless parts of the exact (T"^) on event horizons can reproduced by the auxiliary 
field method based on the trace anomaly. 
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Several approaches to approximating (T*^^) have been discussed previously in the literature, 
developed with important special cases in mind, such as static geometries with a timelike Killing 



field 
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25( 1 ■ These approximations are quite successful for regular states 
in the Schwarzschild geometry, but fail when compared to the numerical results for (T'^) in the 
charged RN spacetimes. Specifically, the previous approximations invariably yield a renormalized 
{T"'f) which grows logarithmically without bound as the horizon of any charged RN black hole is 
approached. For the case of an extreme Reissner-Nordstrom black hole there is an even stronger 
linear divergence. On the other hand direct numerical evaluation of (T'^) in the Hartle-Hawking- 



Israel 



26j thermal state shows no evidence for any of these divergences Q, 3]. We review 



the previous approximation methods and compare them both to the auxiliary field method and the 
direct evaluations of (T'^) in Section 4. 

Our main purpose in this paper is to show that the auxiliary field effective action and stress 
tensor determined by the trace anomaly leads to a practical semi-analytic approximation technique 
which allows for a finite (T'^) on the event horizons of all electrically charged black holes, including 
the extreme Reissner-Nordstrom (ERN) case oi Q = M. Although the stress tensor diverges 
on the horizon for generic solutions of the auxiliary field equations, it is possible to adjust the 
homogeneous solutions of these linear equations to remove the divergences. The ERN case is 
particularly interesting, since its Hawking temperature vanishes and its degenerate horizon structure 
leads to potentially more severely divergent terms in the stress tensor. These leading divergent 
behaviors can be determined analytically by a power series expansion of the auxiliary fields in the 
local vicinity of the horizon, and explicitly cancelled, if desired. 

The paper is organized as follows. In the next section we review the effective action and stress 
tensor of the trace anomaly in the auxiliary field form introduced in [2]. In Section 3, we apply 
the general approximation algorithm to static, spherically symmetric spacetimes, reviewing briefly 
the uncharged Schwarzschild case, and then extending the analysis to the generic charged Q < M 
RN cases, and the ERN Q = M case. We determine in each case the conditions on the series 
expansion coefficients of the auxiliary fields necessary for a regular (T^^) on the horizon. In Section 
4 we solve the regularity conditions and plot the results for the simplest solution in each case, 
comparing them to previous analytic approximation schemes, and direct numerical evaluations of 
(T"^). Section 5 contains our Conclusions, while the Appendix catalogs the complete list of solutions 
to the regularity conditions. 
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II. STRESS TENSOR FROM THE TRACE ANOMALY 

Classical fields satisfying wave equations with zero mass, which are invariant under conformal 
transformations of the spacetime metric, gab e'^'^dab have stress tensors with zero classical trace, 
T"^ = 0. Because the corresponding quantum theory requires an ultraviolet (UV) regulator, classical 
conformal invariance cannot be maintained at the quantum level. The trace of the stress tensor is 
generally non-zero when h ^ 0, and any UV regulator which preserves the covariant conservation 
of T'^, a necessary requirement of any theory respecting general coordinate invariance, yields an 
expectation value of the quantum stress tensor with a non-zero trace: (T'^) ^ 0. This conformal or 
trace anomaly is therefore a general feature of quantum theory in gravitational fields, on the same 



footing as the chiral anoma^ 



y in QCD responsible for the experimentally measured decay of the vr' 







meson into two photons [27[. 

In four spacetime dimensions the trace anomaly takes the general form [2^ , 



(n) = bF + b' (^E-'^nE}j+b"nR + Y, PiHi ■ (2.1) 
In Eq. (12. ip we employ the notation, 

E =*Rabcd = RabcdR"'^'^'^ — ^RabR""^ + R^ i and (2.2a) 

= ^abcd^ — ^abcd^ — '^^abJ^ + ■ (,2.2bj 

with Rabcd the Riemann curvature tensor, *Rabcd = ^^abefR'^^cd dual, and Cabcd the Weyl confor- 
mal tensor. The coefficients b, b' , and b" are dimensionless parameters proportional to h. Additional 
terms denoted by the sum ^ • PiHi in ()2.ip may also appear in the general form of the trace anomaly, 
if the massless field in question couples to additional long range gauge fields. Thus in the case of 
massless fermions coupled to a background gauge field, the invariant H =ti{FahF°'^) appears in 
(j2.ip with a coefficient (3 determined by the beta function of the relevant gauge coupling (29| . 

The form of (j2.ip and coefficients b and b' do not depend on the state in which the expectation 
value of the stress tensor is computed. Instead they are determined only by the number of massless 
fields and their spin via 

h = Y^^^ + + l^A^V-) , (2.3a) 
y = -^^^ (Ns + UNf + 62iVv.) , (2.3b) 
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with Ns the number of spin fields, Np the number of spin ^ Dirac fields, and Ny the number of 
spin 1 fields Henceforth we shall set h = 1, although it should be remembered that any effect 
of the anomaly in which the b, b' and Pi coefficients appear is a one-loop quantum effect. 

The trace anomaly determines the conformal variation of the one-loop effective action of the 
matter fields in a general curved background. A covariant, non-local form of this effective action was 
first given in Ref. Q] . One consequence of the effective action due to the anomaly is that the scalar 
or conformal part of the metric becomes dynamica.1, and its fluctuations provide a mechanism for the 



screening of the cosmological vacuum energy 



30, 



31 



32( 1 . The stress tensor, canonical quantization 



of the conformal degree of freedom and physical states of the quantum conformal factor in the 
absence of the Einstein-Hilbert term were studied in Ref. 33l |. 

That the non-local action of the anomaly could be rendered local by the introduction of scalar 

n r " " 

auxiliary field(s), was noted in Refs. 



effective action were given in 



a a 



],y.Q.a. 

Partial forms of the stress tensor due to this 
l33l |. with the authors of [&] initiating the study of the stress 
tensor obtained from the effective action of the anomaly with auxiliary fields as an approximation 
scheme in the Ricci flat case of Schwarzschild spacetime. The general, complete form of the stress 
tensor in terms of two auxiliary fields in an arbitrary curved spacetime was given in ^. This 
auxiliary field effective action is of the form. 



(2.4) 



with 



in terms of the two scalar auxiliary fields ip and corresponding to the two non-trivial cocycles of 
the Weyl group in four dimensions [7]. 

The effective action ()2.4p - ()2.5p is a spacetime scalar integral over local fields. Hence varying it 
with respect to the metric yields two covariantly conserved stress tensors Eab and Fab, bilinear in 
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the scalar auxiliary fields. Explicitly, these are [4] 

Eab = -2(V(„99)(V;,)nvj) + 2V^[(V,v7)(VaVfeV5)]-^VaVfe[(Vc(^)(VV)] 



and 



+^5af, {-3 + □ [(Ve^)(VV)] + 2 (Sii^"^ - Rg^'') (VeV3)(VdV')} 

2 8 4 

- - Va Vb □ (/P - 4 VcVdV? - 4 Vfe) VcV? +-Rabnip + -RVaVbip 

-l{^{aR)^b)^ + l9ab {2nV + 6i?^^v,Vd(/p-4i?n(^ + (v^fl)v,(^} , (2.6) 

Fab = -2 (V(„V3)(V6) DV) - 2 (V(„^)(Vfc) Dv?) + 2 [(VcV?)(VaVbV') + {Vci^){VaVb^)] 
-^VaVb [(V,(/;)(V^V^)] + ^ (V,V?)(V^V) - 4i?% [(Vfc)(^)(V,V) + (V;,)^)(Ve(/^)] 

+^ii(V(,v9)(Vb)V) + ^5a6{ - 3(nv3)(nV') + □ [(Vc<^)(v^v)] 

+2 (Si?^'^ - (Ve'^)(VrfV')} - 4VeVrf (C(,%V) -2Ca\''RcdV^ 

V(aii) Vb)V+^5afe {2n2^ + 6i?'='^VeVdV-4iZnV+(V'^«)(VcV')} • (2.7) 



2 
3 



These tensors have the local geometrical traces, 



E"^ = 2A4<f = E --DR, (2i 
3 



Fl = 2A4V = F = CabcdC"'" , (2.8b) 

where the latter half of these two equations follow from the independent Euler-Lagrange variation 
of (|2.4p - ()2.5p with respect to the two auxiliary scalar degrees of freedom, if and ^. 



The fourth order scalar differential operator appearing in these expressions is 



A4 = + 2i?'^^VaVfc - ^Rn + ^{V''R)Va = Va (v^V'' + 2R''^ - ^Rg"^^ Vfe . (2.9) 

By solving the fourth order linear equations (12. 8p determined by this A4 for the two auxiliary fields, 
ip and V) ^-nd substituting the results into the stress tensors (j2.6p and ()2.7p we obtain a general 
approximation algorithm for (Tt) for conformal matter fields of any spin in an arbitrary curved 
spacetime. That is, 

(T't) ^ Tt ^P] = + bF^l . (2.10) 
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is an approximation to the exact stress tensor expectation value. Since the dependence of the 
stress tensor Tt[9?, "0] the spin of the underlying quantum matter fields arises purely through 
the numerical coefficients b and b' through (12. 3p . and Pi, there are no new equations to be solved 
for quantum fields of different spin. 

The freedom to add homogeneous solutions of ()2.8p to any given inhomogeneous solution allows 
the tracefree components of the stress tensor (j2.10p to be changed without altering its trace. This 
corresponds to the freedom to change the boundary conditions and the state of the underlying 
quantum field theory without changing its state independent trace anomaly. As shown in U the 
auxiliary fields and traceless terms in the stress tensor (j2.10p generally diverge on event horizons, 
which provides a coordinate invariant meaning to large quantum backreaction effects on horizons. 
These state dependent divergences can be analyzed and removed by specifying boundary conditions 
for the auxiliary fields on the horizon. One then has an approximation scheme for the expectation 
value (Tt) in regular states as well. 

In order to characterize the nature of the approximation (|2.10p . we recall the general decompo- 
sition of the exact quantum effective action into three parts, 



according to its transformation properties under global Weyl transformations [7]. The local action 
Siocai can be expressed purely in terms of local contractions of the Riemann curvature tensor and 
its derivatives. In addition to the classical Einstein-Hilbert action Siocai consists of an infinite series 
of higher dimension curvature invariants multiplied by increasing powers of an inverse mass scale. 
These terms, consistent with a general effective field theory analysis of gravity, give higher order 
geometric contributions to the stress tensor, which remain bounded and small for small curvatures. 
The two remaining terms in ()2.1ip are generally non-local. Any non-local terms involving a non- 
zero fixed mass parameter can be expanded in a series of higher derivative local terms multiplied 
by powers of the inverse mass and regrouped into Siocai- Hence we need consider only those non- 
local terms which are not associated with any fixed mass or length scale in the remaining terms of 
Sexact- These must be either strictly Weyl invariant, denoted by Sinv in (|2.1ip . or break local Weyl 
invariance, yet without introducing any explicit mass or length scale. Up to possible surface terms, 
these are just the geometric terms required by the trace anomaly (|2.ip . The associated terms in 
Sanom scalc logarithmically with distance or energy, and are composed of the two distinct cocycles 
of the Weyl group, given by ()2.4p . 




anom i 




n 
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When the background spacetime is conformally flat or approximately so, the Weyl invariant 
action Sinv may be neglected, since it vanishes in the conformally related flat spacetime, in the usual 
Poincare invariant vacuum state. In that case we expect Sanom to become a good approximation to 
the non-local terms in the exact effective action (j2.1ip . and the corresponding stress tensor (j2.10p 
to become a good approximation to the exact quantum stress tensor, up to well-known local terms. 
Thus the approximation (j2.10p amounts to the neglect of Sinv, or more precisely, those parts of 
Sinv which cannot be expressed in terms of local terms or parameterized by homogeneous solutions 
to the auxiliary field equations (|2.8p . 

Because of the conformal behavior of fields near an event horizon, where the effects of mass 
terms become subdominant, one might expect the leading behavior of the stress tensor (j2.10p also 
to match that of the exact (T^^) in the vicinity of the horizon. In we tested this hypothesis in 
Schwarzschild and de Sitter spacetimes, finding that the freedom to choose homogeneous solutions 
to the auxiliary field equations (j2.8p allows for all possible allowed behaviors of the exact stress 
tensor near the horizon. Indeed for states with diverging stress tensor on the horizon, the anomalous 
stress tensor (I2.10p gives the correct leading and subleading behaviors of the exact (T*^^). 

When attention is restricted to states with regular behavior on the horizon, which from the 
point of view of conformal invariance are subleading with respect to the divergent terms, then the 
stress tensor ()2.10p can often be adjusted to give the exact finite value of (T'^) on the horizon 
as well. However, the global fit of ()2.10p is in only fair quantitative agreement with the numeri- 
cally computed expectation value (T'^) far from the horizon for the Hartle-Hawking-Israel state in 
Schwarzschild spacetime. In regular states the neglected terms in Sexact, which remain bounded 
at the horizon, are comparable in magnitude to subleading terms of Sanom- Hence neglect of Sinv 
is expected to yield a poorer global approximation to the stress tensor of such regular states even 
when it is possible to match the exact behavior of {T'\) at the horizon. 

The Schw..chM and de Sitte. e..e, considered in U are spee.l „. that the anx.lia.y field 
equations (|2.8p may be solved analytically. However, the approximation (j2.10p does not require 
this, and in the following we extend the auxiliary field method to general static, spherically sym- 
metric geometries, focusing specifically on the RN family of charged black holes. These black hole 
spacetimes provide an interesting testbed for the auxiliary field stress tensor, whose qualitative 
and quantitative features may be compared with both previous approximation methods and direct 
numerical evaluation of (T"^). 
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III. REGULAR STRESS TENSORS IN REISSNER-NORDSTROM SPACETIMES 

The effective action and stress tensor of the auxiliary fields (I2.10p is defined in any spacetime, 
regardless of special symmetries. However it becomes particularly useful as a method of approx- 
imating the expectation value of the quantum stress tensor in spacetimes with a high degree of 
symmetry, such as spherical symmetry. The line element for a general static, spherically symmetric 
spacetime can be expressed in terms of two functions of the radius in the form, 

ds" = -f{r)dt' + ^+r''dn\ (3.1) 
n{r) 

Assuming that the state in which we evaluate the the stress tensor is also spherically symmetric 
and stationary in time, we may make a static, spherically symmetric ansatz for the auxiliary fields 
as well, i.e. 

(p = if{r) (3.2a) 
ip = V(r), (3.2b) 



so that the equations (|2.8p become ordinary differential equations in r. In some cases it is possible 
also to add terms with linear time dependence to the auxiliary fields, i.e. (p = ip{r) + at and and 
ip = ip(^r) + a't in order to allow for the possibility of non-vanishing T''^ flux components which are 
also independent of the Killing time t. Higher powers of t or more complicated time dependence in 
the auxiliary fields lead to non-stationary stress tensors. 

For a general static, spherically symmetric spacetime, with the fields in a spherically symmetric 
quantum state, the stress tensor is given by its three independent diagonal components, 

T\ = -p{r) (3.3) 
r; = p{r) (3.4) 
T% = p^{r), (3.5) 

together with a possible non-zero off-diagonal fiux component T^^. These components obey the 
covariant conservation conditions, 

Van = ^ + :^fiP + p) + -{p-P±)=0 (3.6) 
ar I J ar r 



and 



v.n = -L-^(V^r';) = o. (3.7) 
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Eq. (j3.7p can be trivially integrated to obtain 

with the integration constant L the luminosity of a localized source. 

The approximation of (|2.1U|) can be applied to arbitrary spacetimes, spherically symmetric or 
not, with or without an horizon. For definiteness we restrict our attention in this paper to the 
family of Reissner-Nordstrom spacetimes with the equal metric functions, 

/(r) = h{r) = 1-^ + 91^ (3.9) 
depending upon M, the mass and the electric charge. We then can distinguish three cases: 

(i) . Q = 0, Schwarzschild spacetime; 

(ii) . < Q < M, generic Reissner-Nordstrom (RN) spacetime; 

(iii) . Q = M, extreme Reissner-Nordstrom (ERN) spacetime. 

We shall discuss each of these cases in detail separately. 

Since (|3.9p is quadratic in 1/r there are two values of r at which /(r) vanishes. These are 



r± = M± - Q2 _ (310) 

When Q = 0, r+ = 2M is the usual Schwarzschild horizon. In the ERN case of Q = M, the 
two values r± coincide and the horizon becomes degenerate. The character of the spacetime, the 
solutions to the auxiliary field equations and the corresponding stress tensors derived from them 
are quite different in each of the three cases. Henceforth we reserve the designation RN for the 
generic, charged black hole solution of case (ii). 

For conformal field theories, the trace of the stress tensor is given purely by the trace anomaly. 
This provides us with another relation for the diagonal components of the tensor, namely 

- p + p + 2p± = T = b'E + bF (3.11) 

Furthermore, defining 

Q{r)^P±-'^ (3.12) 
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and integrating (j3.6p with (j3.1ip gives [3j] 

Pi-) = M - rV) e^^jy-<giT^y^, (3.13) 

where C — L denotes an overall constant of integration. The only unknown function in this 
expression is 0(r). Thus it is sufficient to compute 0(r) from the auxiliary field stress tensor to 
determine all the diagonal components of the stress tensor, p{r), p±{r) and p{r) from equations 
(I3T3I1 . (I3T2I1 and (ISTTTl respectively. From Eq. i^B) with h{r) = f{r), the off diagonal flux 
component is 

Since the curvature invariants in the trace (j3.1ip remain finite as the horizon at r = r-|_ is 
approached, it is clear from (j3.12p that p± is finite on the horizon, provided 0(r) remains finite 
there. Moreover if /(r) possesses an isolated simple zero at r = r_|_, vanishing linearly in 

s = (r — r+)/r+. (3.15) 

as r — > r-i-, so that the event horizon is non-degenerate, then (j3.13p shows that a divergent term 
in the other components of the stress tensor can arise only if C — L ^ 0. Hence 0(r) remaining 
finite as r — > r_|_, and C = L, are necessary and sufficient conditions for finiteness of all components 
of {T"^) in the static Killing frame of (j3.ip for the non-degenerate Reissner-Nordstrom horizons 
(0 < Q < M). 



Finiteness of the stress tensor in the frame of a freely falling observer requires also that 

IT I 1 



P 4/2 



{P + P)f + J^ <oo- (3.16) 



27rr2 

If 0(r) is regular at r = r+ and can be expanded in a Taylor series in s near s = 0, then it is easy 
to see that the condition (j3.16p is satisfied automatically if C = 0. In the case of the past horizon 
one should demand finiteness of \Tyy\/ p instead of (j3.16p . This forces the integration constant to 
be C = 2L instead. Thus regularity on both the past and future horizon requires Q{r^) < oo and 
C = L = 0. In this case the regularity condition (j3.16p is equivalent to the condition 

= y|p + p| < oo . (3.17) 

In some cases, i.e. those with logarithmic terms in the auxiliary fields, 0(r) is not analytic at 
r = r_(_ and cannot be expanded in a Taylor series there. In those cases, the finiteness condition 
(|3.16p will give additional conditions on the behavior of the auxiliary fields on the horizon. 
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Before restricting our attention to the static, completely regular class of stress tensors, let us 
emphasize that this is not the generic case, as a different choice of the free integration constants C 
and L would lead to different physical behavior on the horizon. It is well known that the coordinate 
singularity of the metric (j3.1|) at the horizon r = r+ where f = h vanishes may be removed by 
a (singular) coordinate transformation, producing the complex analytic extension of the Reissner- 
Nordstrom geometry [35]. Although any regular transformation of coordinates is allowed by the 
Equivalence Principle, and cannot lead to physical effects, sin^it/ar coordinate transformations, like 
singular gauge transformations in gauge theory, must be treated with some care. New topological 
configurations such as monopoles or vortices are associated with such singular gauge transforma- 
tions. Thus, although the complex analytic extension of a black hole spacetime may seem quite 
natural mathematically, analytic continuation actually involves a physical assumption, namely that 
there are no stress tensor sources to the Einstein equations localized on or near the horizon. Because 
of the hyperbolic nature of Einstein's equations such stress sources with effects transmitted along a 
null surface are perfectly allowable, even classically. When the expectation value of the stress tensor 
of quantum fields is considered, with its sensitivity to the wavelike, non-local coherence effects of 
quantum matter, the assumption of analyticity on the horizon is not at all automatic, and is not 
required by any general principles of quantum theory. 

The effective action of the conformal anomaly, and its associated auxiliary fields indicate that 



non-regularity of the stress tensor on the horizon is to be expected in the generic case as well [2]. 
A tuning of the integration constants of the solutions of the linear equations (j2.8p for ip and ip is 
necessary to prevent 0(r) from diverging as r — > r_(_, with the generic behavior of all the stress 
tensor components near the horizon being proportional to as / — > in the Schwarzschild case. 
Since the auxiliary fields are spacetime scalars, it is clear that this behavior is in no contradiction 
to the Equivalence Principle. In fact, the divergences have a perfectly coordinate invariant origin 
in terms of the homogeneous solutions to the equations ()2.8p 



where i^T = ^ is the Killing field of the static Schwarzschild or Reissner-Nordstrom geometries, 
timelike for r > r^. This defines the rest frame of the configuration, which is independent of 
coordinate redefinitions. The allowed divergences in the stress tensor on the horizon as /(r) 
and ip^tlj ^ zbco are related therefore to the behavior of this Killing invariant of the global geometry 
becoming null. The state of the quantum matter fields, specified on a complete Cauchy surface of 



n 




(3.18) 
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the spacetime is necessarily defined in a non-local way, and hence expectation values of T*^ can be 
sensitive to the divergences of (|3.18p at this null surface, notwithstanding the finiteness of the local 
Riemann curvature at r = r_|_. 

By choosing suitably restricted solutions of the auxiliary field equations (j2.8p it is possible to 
cancel the In / behavior in (/? and thereby guaranteeing that 0(r-|_) is finite. When the electric 
charge satisfies < Q < M, i.e. excluding the uncharged Schwarzschild and maximally charged 
ERN cases, the stress tensor will have subleading In s and In^ s divergences as s — > as well. Then 
@{r) is not expandable as a Taylor series around r = r+ and the logarithmic terms do not drop 
out entirely from the condition (j3.16p . This leads to additional conditions on the coefficients of the 
expansion around s = for a fully regular stress tensor. 



A. Schwarzschild Spacetime 



In the Schwarzschild case, /(r) = 1 — 2M/r, and the fourth order linear equations 
integrated explicitly for auxiliary fields which are functions only of r. The result is 
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r{r-2M) 2M \2M 

in terms of the three dimensionless constants of integration, c^, Cqo, and q. A fourth integration 
constant would be introduced by integrating (j3.19p once further, but as the stress tensor in the 
Schwarzschild case depends only upon derivatives of ip, a constant shift in (/? plays no role in this 
case. 

The role of the three integration constants appearing in (j3.19p is best exposed by examining the 
limits, 
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oo. (3.20b) 



4M2 4M r 3r2""""V2My 
Hence controls the leading behavior as r approaches the horizon, corresponding to the homoge- 
neous solution of ()3.18p . It is this leading behavior that gives rise to the generic behavior of 
the stress tensor as r — > r_|_. The second integration constant Cqo controls the leading behavior of 
ip{r) as r — > oo, which is the same as in flat space. Non-zero values of Cqo correspond to non-trivial 
boundary conditions at some large but finite volume, such as may be appropriate in the Casimir 
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effect, or if the black hole is enclosed in a box. The constant q is the topological charge of the auxil- 
iary field configuration, associated with the conserved current generated by the Noether symmetry 
of the effective action (12. 4p . ijj ^ ijj + const. Qj. It is responsible for the Inr terms in (I3.20bp and 
the corresponding stress tensor (j2.6p . 

To the general spherically symmetric static solution (13.190 we may add also a term linear in t, 
i.e. we may replace ^{r) by 

(^(r,t) = (^(r) + ^t, (3.21) 

with rj an additional free constant of integration. Linear time dependence in the auxiliary fields 
is the only allowed time dependence that leads to a time-independent stress-energy, and this only 
in the Ricci flat Schwarzschild case. The solution for -0 = 'ip{r,t) is of the same form as (j3.19p 
and (j3.2ip with four new integration constants, d^,doo, q' and rj' replacing c^,Coo, q, and rj in 
ip{r,t). Adding terms with any higher powers of t or more complicated t dependence produces a 
time dependent stress-energy tensor. The stress tensor (I2.10j) does not depend on either a constant 
ipo or V'o- 

The stress-energy diverges on the horizon in an entire family of states for generic values of the 
eight auxiliary field parameters (c^ , q, Coo,r]; , q' , doo,r]'). Hence in the general allowed parameter 
space of spherically symmetric macroscopic states, horizon divergences of the stress-energy are quite 



generic, and not restricted to the Boulware state [36]. In addition to the leading behavior, there 
are subleading s^^, In^ s and Ins divergences in general. It turns out that only three of these four 
are independent, and the three conditions. 



-{b'cH + 2bdH)cH + il{b'ii + 26??') = [s-"^ in (3.22a) 

{q - 2)[b'{q - 2) + 2b{q' - 2)] = (in^ s in (3.22b) 

b[{q - 2)(18doo - -i^H - 40) + [q' - 2)(18coo - 30ch - 40)] + 

+b'{q - 2) (IScoo - 30ciT - 40) = (in s in (3.22c) 

are all that are required to remove the divergent behaviors in 0, indicated in parentheses as s ^ 0, 
including a possible divergence at the horizon. The first of these conditions, (|3.22ap corrects 
a sign error in Eq. (5.14b) of Ref. [2I (where the notations p,p' were used in place of the present 
r/,r/' for the parameters of the linear time dependent terms in -0 respectively). 
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The regularity condition p.l6p on T'^'^ gives in addition, 

b (^^^ + 2cH + 2^^) + b' + 2^^) = b'm + b{q7]' + q'r]) = ^ (3.23a) 

in T"") 

6[ci^ (g' - 2) + dj^ - 2)] + t'dniq - 2) = (In s in T"") (3.23b) 
If we are interested in strictly static, regular states with C = L = in (|3.13p and ()3.14p . then we 
obtain an additional condition, namely the quantity in (|3.23ap proportional to the luminosity must 
vanish. The condition ()3.23bp eliminates a possible subleading logarithmic divergence in Tuu/ f'^ 
on the horizon. It is clear that the choice q = q' = 2 satisfies this condition as well as ()3.22bp and 
(j3.22cp . This illustrates the general property that when logarithmic terms are taken to vanish on 
the horizon, the conditions that 0(r_|_) be finite on the horizon and C = L = are sufficient to yield 
a fully regular stress tensor on the horizon, and the number of independent conditions is reduced. 



B. Generic Charged RN Spacetimes 

The generic Q > charged Reissner-Nordstrom spacetimes are not Ricci flat and an analytic 
solution of the fourth order equations (j2.8p in closed form no longer appears possible. In addition, 
another difference from the Schwarzschild case, stemming from the same non- vanishing of the Ricci 
tensor when Q > 0, is that the linear time dependence (|3.2ip in the auxiliary field ip now produces 
time dependence in the stress tensor (|2.7p . and hence is disallowed for a static state. Linear time 
dependence in ip is allowed, leading however to a non-zero fiux. Hence for strictly static states we 
must set r] = r]' = 0. However, the non-Ricci fiatness means also that a possible constant term 
in if which drops out of the stress tensor (j2.10p when R^t = now survives as a non-trivial free 
parameter in this case. Hence we have seven remaining integration constants in all for spherically 
symmetric static auxiliary fields (j3.2p in the charged RN case. 

Since the fourth order differential operator, A4 involves a total derivative, cf. Eq. (j2.9p . equa- 
tions (j2.8p can be integrated once, to obtain the second order equation for (p'{r), 

-2Q2j^(^' = eo + ^ / r'^drE, (3.24) 

where cq is an integration constant. Here we have used the facts that 

R\ = -^ (3.25) 

and = in a general RN spacetime. 



dr 



dr^'^ ^ 
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A power series representation of the general solution of the second order equation (j3.24p for ip' 
in powers of s = (r — r_|_ ) /r_|_ is easily derived. Since the right side of (j3.24p is regular at s = 0, it 
can be expressed as a Taylor series in the form, 

eo + ^ r r'drE = Y,enS^ . (3.26) 

^ •^^+ n=0 

To find the leading behavior of near s = let ~ s"^ for some 7, and thereby obtain from (j3.24p 
and (|3:26]) . 



£27(7 + 1)5^ + 2 £7(7 + 2)^5^+1 



■)2 



+ U7 + l)(7 + 2)-^ + 2£ 



-(7 + 1)' + 



2Q^ 



I + 0(^7+3) 



5]e„s", (3.27) 



n=0 



where 



2v/M2 - g2 



(3.28) 



r+ M + v^M2 - g2 

in the general case. From p.27p we observe that 7 = — 1 gives the most singular behavior allowed 
for 99' as s ^ for general eo 7^ 0, and > r_. Note that 99' ~ s"^ agrees also with the leading 
behavior in (j3.20ap for the exact solution in the uncharged Schwarzschild case. 

Because the singular behavior is allowed for (/?' for all < Q < M, (/? has at most logarith- 
mically singular behavior near the RN horizon, and its general series expansion is of the form, 

^{r) = (^ns"" + ^ns"" In s . (3.29) 

ri,=0 n=0 

Note that the general logarithmic behavior near s = is that expected from (j3.18p on geometrical 
grounds, since e~^''- with iph = Ins is the conformal transformation needed to bring the RN metric 
to the ultrastatic optical metric, and remove the singularity caused by the Killing field changing 
character from timelike to null. 

Substitution of the expansion ()3.29p into the differential equation (I3.24p gives recursion relations 
for the coefficients {an,^n}- The set of four coefficients, (oq, ai,£o;^i) are free integration constants 
for (p, parameterizing the general solution of the fourth order equation, with all higher order an>i 
and 6n>i determined by the recursion relations in terms of these four integration constants. In the 
Schwarzschild limit, £0 ch,(i —>■ q — 2,ai — >Ciy + g+ | — 3coo respectively. The constant cq is 
also determined in terms of £0 by (j3.24p or (j3.27p to be 

02 

eo = -2£4-V. (3.30) 

r± 
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In like manner substitution of the series expansion, 




(3.31) 



n=0 



into the equation for tp shows that (cq, ci, Aq, Ai) are the four free integration constants parame- 
terizing the general solution of (|2.8b|) . However, since a constant shift in ip is an exact symmetry 
of the action (j2.4p . cq does not appear in the stress tensor (|2.10p . and we are left with the seven 
effective free integration constants, (ao, ai, ^O) ^i! ci, Aq, Ai) in total. 

The general form of the divergences of the anomalous stress tensor at the RN event horizon 
may be found from the power series expansions of the auxiliary fields there. In the T^g component 
there are s~^, and In^ s divergences. The leading behavior in all components of T'^ in 
coordinates (13. ip may be understood from the approximate conformal symmetry which applies near 
the horizon and the conformal weight of the stress tensor. Namely, since the near horizon geometry 
is conformally flat, with e~'^'*^*/''+ the conformal transformation to locally flat space, distances 
scale like e"^" = /a, while energies scale like f~2, and energy densities like ~ near the 
non-degenerate RN event horizon. 

Requiring the coefficients of the leading and subleading In^ s and In s possible divergences in 
Q to be zero gives three conditions on the integration constants of the auxiliary fields. A possible 
divergence in Q turns out to be linearly dependent on the first three, and is canceled automatically 
when these three conditions are imposed. However from (j3.13p with L = there remains a possible 
divergence in T^g unless C = 0. This is not automatic and gives a fourth condition. A fifth 
and final condition comes from the necessity of canceling the In s divergence in the freely falling 
frame, T**" of (j3.17p . A possible In^ s term in T"" drops out after we have satisfied the ffi'st three 
conditions, and gives no further condition. Hence we end up with five algebraic relations for the 
seven constants of integration, viz., 



it'll + 6Ai)[3 - e(ai + ^i) + 2(3e - 1)^] + 6£i[3e - e(ci + Ai) + 2(3e - l)Ao] =0 (Ins in T^^g) (3.32c) 
b{6eio + 6Ao + 4e^iAo + 4e4Ai - eAi) + 6'(64 - eh + 4e4^i) = (C = 0) (3.32d) 
6[18e(l - e)io + 9e(e - + 3(1 - 4e + 3e^)Ai - (3 - 12e + 20e'^){iiXo + ^Ai)] 



where e is given by (I3.28p . The solutions of these conditions will be discussed in the next section. 




(3.32b) 



(3.32a) 



+6'[3(1 - 4e + 3e^)£i - (3 - 12e + 20e2)£i4] = (In s in T"") 



(3.32e) 
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C. ERN Spacetime 

When Q = M, = = M, the horizon becomes degenerate, and the RN metric function 
/(r) goes to zero quadraticahy as r ^ r-t, (s ^ 0). This quite different behavior of the spacetime 
near the horizon is reflected also in the behavior of the conformal differential operator A4 and its 
solutions. Referring back to (j3.27p we observe that when e = the first two terms vanish identically, 
and the coefficient of the s'^'^'^ term becomes the leading one, with coefficient (7 + 1)(7 + 2) — 2 = 
7(7 + 3). Thus the more singular behavior 7 = —3 is allowed by (|3.27p . and the structure of the 
divergent terms in the solutions to the auxiliary field equations and stress tensor (j2.10p becomes 
quite different in the ERN limit. For this reason the number of conditions necessary for the stress 
tensor (j2.10p to remain finite increases, and the conditions become more stringent. 

Because of the leading behavior allowed by (j3.24p for 99' and tp', the power series expansions 
of the auxiliary fields in the ERN case are of the form, 

ip = ^ ans"" + ^ ens'" Ins (3.33a) 

n=-2 n=0 

^ = Cns'^ + ^A^s^lns (3.33b) 

n=-2 n=0 

instead of ()3.29p and ()3.3ip . Substitution of these series into the differential equations for ip and -0 
shows that (a_2, a_i, ao, oi) and (c_2, c_i, cq, ci) are eight free integration constants, with cq again 
playing no role. The logarithmic terms and all the other coefficients are related to this set of seven 
coefficients by recursion relations. Since the leading behaviors of p and tp are now as s ^ 0, 
the stress tensor (I2.10p gives T^g with s~^,s~^, s^^, s"^, s~^lns, In^ s and Ins divergent terms. 
The leading behavior again can be understood from the conformal weight of the stress tensor 
under conformal transformations near the horizon, since /^^ ~ in the ERN case. 

It turns out that the conditions for removing the leading and all subleading divergent behaviors in 
the auxiliary field T'\ are not independent, and only /oitr independent conditions on the integration 
constants are sufficient. We can choose these four to be: 

a_2(26c_2 + b'a-2) = (s~^ in 9) (3.34a) 

6(-10a_2C_2 + a-2C-i + a_ic_2) + 6'(-5a?_2 + a_2a-i) = (s"^ in 6) (3.34b) 

b [3c_i - 2a_2(9 + 3ci + 4c_i - 64c_2) - c_2(36 + 6ai + 8a„i)] + (3.34c) 

+b' [3a__i - a_2(36 + 8a_i + 6ai - 64a_2)] = (s^^ in 9) 

a_i(25c_i + b'a-i) = (Ins in 9) (3.34d) 
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In particular the last of these conditions removes all logarithmically divergent terms from Q. 
The requirement (13.170 that T"" be finite gives two additional conditions, viz., 

b{-15 + 15a_i - 96a_2 + 13c_i - 104c_2) + 

+&'(13a„i - 104a_2) = (s"^ in T™) (3.35a) 
b [a_2(486 - 32c_i) + a_i(3ci - 72 + 8c_i - 32c_2) + 27 + 3c_i(ai - 10) + 288c_2] 

+b' [Aaii + 288a_2 - a_i(30 - 3ai + 32a_2)] =0 (Ins in T"") (3.35b) 

Hence we have six independent conditions on the seven active integration constants of the auxiliary 
fields in order to obtain a fully finite stress tensor from the trace anomaly in freely falling coordinates 
in the ERN case. Despite the discontinuously singular behavior of the geometry as Q ^ M, the 
more singular behavior of the auxiliary fields in this limit, and the greater restrictiveness of the 
finiteness conditions (j3.34p - (j3.35p compared to (j3.32p . it is still possible for the approximation 
algorithm for the stress tensor based on the anomalous effective action (j2.4p to yield a fully finite 
stress tensor on the ERN horizon. This is qualitatively different from all previous approximation 
schemes. 



IV. COMPARISON WITH PREVIOUS APPROXIMATIONS AND NUMERICAL 

RESULTS 



Since the method of computing the stress tensor from the anomaly action with the use of 
auxiliary fields is relatively new, it is interesting to compare it to previous approximation methods, 
as well as direct numerical evaluations of the renormalized {T'^f) whenever the latter are available. 
An analytic approximation to (T'^) for thermal states of conformal fields in non-conformally flat 



static spacetimes was derived by Page 



19l|. Brown and Ottewill [2^, and Brown, Ottewill and 



Page [21[. T 
heat kernel 



lis approximation is based upon the semi-classical approximation to the proper time 



37( 1 , and the properties of the exact one-loop effective action under the conformal 



transformation , 



2uj ~ 
gab = e Qa 



(4.1) 



For a classical conformal field the dependence of the exact effective action upon lo is determined 
completely by the trace anomaly in the form, 



Sexactis] — ^1 



exact I 



+ bA[Lo;g] + b'B[u;;g] 



(4.2) 
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where A[uj; g] and B[uj; g] are known (non-linear) functionals of uj and gab which are given in 20|,|21|. 
If the conformal transformation e^^ is chosen to be f{r), for the static, spherically symmetric line 
element (I3.ip . then the conformally transformed metric gab becomes the ultrastatic, optical metric, 
for which gu = —1- If, in addition, the original physical metric gab is that of a static Einstein space 
(i.e. one for which Rab = ^Sab), then the invariants appearing in the anomaly. 



-g 



■■ dab ■ 



5A 

Sgab 
5B 



0, 



9=g 



E 
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V -9' Saab "J=3 

vanish for the conformally related ultrastatic metric c/ab, and the A and B terms in 



(4.3a) 

(4.3b) 

reproduce 



the correct trace for the physical metric gab- Thus, with this choice of lv it is consistent simply 
to neglect Sexact[g\ in (|4.2p . and approximate the full (T'^j,) by the terms coming from yl[a;;g] and 
B[oj;g] which are known analytically. 

Because it is also based on the form of the trace anomaly, the Page-Brown-Ottewill (PBO) 
approximation is related to the approximation scheme based on the auxiliary field effective action 
of the anomaly. Indeed, by making use of the conformal transformation property, 

2 



-a 



R + Q{Uu - g^^daUjdbOJ 



(4.4) 



it is not difficult to show that 

hA[uj]g] + h'B[u}]g\ = -Twz[g; -w] + 

where 



+ h') , 4 

a x{y—gR 
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-gR^ 



(4.5) 



^Wz[g] -Uj] = Sanom[g\ " Sanoni[g = e ^'^fi'] (4.6) 

is the Wess-Zumino effective action for the anomaly obtained in P, 7|. Since T^z is a quadratic 
functional of uj, the relation (j4.5p shows that the complicated cubic and quartic terms of the PBO 
effective action are simply the result of adding an R? term to T\yz, with a corresponding b"\I\R 
term in the trace anomaly (j2.ip . The PBO effective action Vl^pBob] is related then to the anomaly 
effective action Sanom of (j2.4p by 

+ 6') 



WpBo[g\ = WpBO 
~^ Sanom [g] 



+ 18 

ib + b') r,4 



~gR^ 



Sn 
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(4.7) 
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This consistent with the decomposition (j2.1ip . in which all terms that depend only upon g^b are 
viewed as invariant under the local Weyl transformation ()4.ip . Because of the different linear com- 
bination of invariants in (|4.3p from the b and b' terms in (12. ip . the PBO conformal transformation 
is given by an u; which satisfies 

iA^LO = E - F --nR, (4.8) 
3 

and therefore corresponds to a particular linear combination (i.e. (p — ip) of the auxiliary fields, 
with a particular choice of homogeneous solution. 

Although the PBO approximation is related to the one based on Sanom, there are two important 
differences. First, in the PBO approximation the conformal transformation e'^'^ to the ultrastatic 
metric is fixed up to a linearly time dependent term in uj, by the requirement that the invariants in 
(j4.3p vanish in the conformally transformed spacetime. Hence there is only one parameter, namely 
the coefficient of linear time dependence free in uj, which is far fewer than the seven integration 
constants corresponding to the freedom to adjust the state dependent and Weyl invariant part 
of the effective action in the auxiliary field method. Secondly, and even more importantly, the 
invariants in (|4.3p vanish in the conformally transformed spacetime only in rather special cases, 
such as if the original metric gab is that of a static Einstein space. Although Zannias suggested a 



modification of the PBO ansatz to account for the correct trace in non Einstein spaces [2^, the 
non- vanishing of the trace in the conformally related space with metric gab means that the original 
PBO rationale for ignoring 5'e^act[5'] (14. 2p no longer applies. Hence there is no especially good 
reason why this modification of the PBO approximation should be accurate, or even finite on the 
event horizon. If one tries instead to find a spacetime where the conditions ()4.3p for the vanishing 
of the trace are satisfied, then one is faced with solving two non-linear conditions for a single to. 
Hence one would not expect that a simultaneous solution of both conditions (j4.3p exists at all for 

'ab- 

In the auxiliary field method, these difficulties are removed completely. The two independent 
ip and "0 fields satisfy the linear equations (j2.8p . for which solutions are guaranteed to exist and 
can be found in any spacetime. Although e~'^ may be regarded as a conformal transformation to 
a spacetime where E — ^\ZIR vanishes, there is no such conformal interpretation for the second 
auxiliary field ip, and none is required. The auxiliary scalars and their stress tensor simply encode 
the same information about the full non-local trace anomaly in a local, generally covariant form, 
and there is no conformally related gab enjoying a privileged status over any other. This is clear 



gab = e '^'^gab, for general g, 
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also from the freedom to add arbitrary homogeneous solutions to ()2.8p . corresponding to different 
choices for the conformal image Qab, and different Weyl invariant parts of the effective action. This 
additional freedom in the linear system ()2.8p makes an exploration of a much wider class of states 
possible with the auxiliary field algorithm, with different state-dependent traceless contributions 
to the stress tensor, all consistent with the correct trace anomaly. It is of course the much larger 
parameter space available in the auxiliary field method which makes it possible to find regular 
stress tensors on both the RN and ERN event horizons. 

The PBO approximation was later rederived in a different way by Frolov and Zel'nikov (FZ), 
by carrying out an analysis of possible terms in the effective action in general spacetimes with a 
static Killing vector field, Ka — — 



at 



231 ]. Although rather different in methodology, the FZ approach 
also applies only to static spacetimes and makes use of the conformal transformation properties 
(conformal weights) of the various terms in the stress tensor. As a result, although the FZ approach 
is not limited a priori to Einstein spaces, the effective action they obtain is in fact equivalent to that 
of PBO (as modified by Zannias), up to local F = Cabcd.C^^'^'^ and the {^J^-g^f' - ^/^R^) terms 
appearing in (14. 5p , which are allowed to have arbitrary coefficients unrelated to the anomaly. These 
particular Weyl invariant terms are mildly behaved on the event horizon, and hence cannot be used 
to cancel any divergences present in the PBO approach. Hence when applied to the non-Ricci flat 
Reissner-Nordstrom geometry, the FZ approximation suffers from the same limitation as that of 
PBO, namely, both predict a logarithmically divergent stress tensor on the RN event horizon, and 
both a linear and a logarithmic divergence on the ERN horizon. 

From our present vantage point, the interesting feature of the FZ approach is the central role 
of the static Killing field Ka- The FZ approach underlines the fact that (T'^) generally depends 
upon global invariant functions of Ka as in (j3.18p . in addition to strictly local invariants such as 
F, E and E?. Since global invariants such as (|3.18p may diverge on the event horizon, and (T'^) 
generally is a function of these invariants, explicitly so in the FZ approach, it is clear that requiring 
such divergences to be absent in (T'^) is a dynamical restriction on the quantum state, not at all 
required by general coordinate invariance. 

It was shown by Howard and Candelas [l^, [3] for the conformally invariant spin field in 
Schwarzschild spacetime that the WKB approximation can be used to write the stress-energy tensor 
in terms of the PBO approximation for each field plus a term containing mode sums that must be 
computed numerically. A similar result was obtained for the massless spin 1 field by Jensen and 
Ottewill [14] . Anderson, Hiscock, and Samuel [l^ (AHS) showed that if the WKB approximation 
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in the high frequency hmit is used for the radial modes of the Euchdean Green's function for a 
massless scalar field with arbitrary curvature coupling in a general static spherically symmetric 
spacetime, then a conserved stress-energy tensor results which in the case of conformal coupling 
has a trace equal to the trace anomaly. For the case of conformal coupling this stress-energy tensor 
is equivalent to that derived by FZ if the three arbitrary constants in their derivation are set to 
zero. A similar approximation was derived by Groves, Anderson, and Carlson 2^ (GAG) for the 



massless spin 1/2 field. In this case the approximation is equivalent to that of FZ if their arbitrary 
constants have the values qf^^ = q^^ = and q^"^ = 1/144. 

nation of the stress tensor for a conformally invariant scalar field in the PBO/FZ 



Huang's eva 

approximation [3^, apparently correcting an error in 22] for the ERN case, shows this same loga- 
rithmic divergence of (j3.16p for Q < M, which becomes a combination of a linear and logarithmic 
divergence in the ERN case.^ Not surprisingly the same divergences were found for the AHS and 
GAG approximations [3]. Thus all pre-existing analytic approximations to (T'^) diverge at 
least logarithmically in the general RN spacetime, and linearly in the ERN case. 

A. Schwarzschild Spacetime 

The six conditions (I3.22|) and (13.230 guaranteeing finiteness of the stress tensor on the 
Schwarzschild horizon can be satisfied in several different ways. The simplest possibility is [3] 



{2b + b')cl +r]{2br]' + b'r]) = 0, (4.9a) 
bd^ = -{b + b')c^ , (4.9b) 
q = q' = 2, (4.9c) 
b'rjq + b{qr]' + q'r]) = , (4.9d) 
which requires fixing only five parameters to satisfy all six conditions (j3.22p and (j3.23p . This is 
because the choice q = q' = 2 eliminates all the logarithmic terms in the auxiliary fields at the 
horizon (hence all In^ s and Ins terms in T*^) and simplifies the remaining conditions considerably. 
The three parameter subset of the original eight parameter family of spherically symmetric auxiliary 
fields is certainly not generic, but still leaves considerable freedom to fit the finite values of the 



^ The third footnote of Ref. [s^ gives an incorrect criterion for finiteness of the stress tensor on the horizon in freely 
falling coordinates, replacing one factor of /(r) in (|3.16p by f^{r). For this reason the conclusions drawn in 
are not warranted by the evaluation of T% given. 
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stress tensor at r = 2M and/or r = oo in the regular Hartle-Hawking- Israel state 



2g]. This is the 

choice which was studied in some detail in The other possible ways of solving the finiteness 



39(] state may be 



conditions are listed in the Appendix. Approximate stress tensors for the Unruh 
obtained by replacing the L = condition (j3.23ap by its correctly normalized value. 

The PBO/FZ approximation works quite well for the conformally invariant scalar field in the 
Hartle-Hawking-Israel state in Schwarzschild spacetime, although the method based on the auxiliary 
field effective action of the anomaly accounts for the behavior of the stress tensor more accurately 
in states with divergent stress tensors on the horizon, such as the Boulware state in Schwarzschild 
spacetime. For the spin 1/2 and spin 1 fields the PBO/FZ approximation works less well [l^. [l8|. 

Balbinot, Fabbri, and Shapiro (BFS) were the first to use two auxiliary fields to study 
the properties of the approximate stress-energy tensor based on the anomaly in Schwarzschild 
spacetime. Thus their approach comes the closest of any previous one to the present work. Their 
two auxiliary fields (0, ip) are related to ours by the linear combinations, —\^—b'ip + btp/^f—b' and 
6/ ^J—h'il) respectively. However, BFS applied conditions to each of these two linear combinations 
of auxiliary fields separately, rather than searching for the more general solutions of the finiteness 
conditions (j3.22p and (j3.23p on the Schwarzschild horizon. This amounts to fixing the traceless 
terms in the stress tensor with fewer free parameters than in the present approach. For this reason 
they were not able to reproduce some features of the exact stress tensor in the Hartle-Hawking-Israel 
state, such as the correct value of (T'^) on the Schwarzschild event horizon or at infinity. In the 
case of the present formulation, after all four finiteness conditions (|3.22p - (j3.23p are satisfied, there 
remains a three parameter family of finite stress tensors. Hence both the values on the horizon and 
at infinity can be adjusted to their correct values by a suitable choice of the integration constants. 

As pointed out by BFS, the subdominant terms in the stress tensor at infinity are also not 
reproduced by the stress tensor of the anomaly, and this undesirable feature persists in our approach. 
Clearly this is because the anomaly action is not equal to the full quantum effective action, differing 
from it by Weyl invariant terms, as in (j2.1ip . In regular states these give rise to additional traceless 
terms in the stress tensor, which would be expected to be of the same order as those in (j2.10p . One 
possibility for improvement is to add the Weyl invariant term, 

Sinv[9;x] = k j ^(fx {-]pc^^X + F:^ (4.10) 

to the total effective action. If we added this term to Sanom, we would have a third auxiliary field, 
denoted here by with four more integration constants to serve as free parameters in the stress 
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tensors of the RN and ERN spacetimes. Although this may allow for more accurate fitting of the 
numerical results for {T'^f) in the RN and ERN cases, we do not pursue this possibility here, and 
thus set /c = in our approximation. 



B. General RN Spacetime 

When the PBO/FZ, AHS, and GAG approximations are applied to the Q > RN spacetimes, 
they predict that (T^^) diverges logarithmically on the RN event horizon for Q < M and linearly and 
logarithmically for Q = M. The reason for the divergence of the PBO/FZ approximation is that 
the semi-classical approximation for the heat kernel of [s^] fails as r — > r+. Thus for massless fields 
all previous approximations are both highly specialized to certain specific classes of spacetimes, 
and also lead inevitably to divergences on the event horizons of non-Einstein static spaces, such 
as the Reissner-Nordstrom metric. The present approach based on the auxiliary field form of the 
effective action of the anomaly does not rely on a WKB approximation, and hence does not suffer 
from these limitations. 

The direct method of evaluating (T'^) from the wave equations of the underlying conformal 



field theory has been carried out for both scalar and spinor fields in RN spacetimes with < 

nnn 

Q < M jl6|, 117|, [18!] • These direct evaluations show no evidence of a divergence as the horizon is 
approached. In the following plots we show the results of these exact numerical evaluations of (T^^) 
and the approximation of the present paper for an RN spacetime with \Q\/M = 0.99. 

The approximation based on the auxiliary field effective action is qualitatively better than that 
of PBO/FZ in all Q > RN regular states in that it allows for a finite stress tensor on the event 
horizon, showing roughly comparable features to those already encountered in the Schwarzschild 
case in Qj. The main difference from the Q = Schwarzschild case to the < Q < M RN cases is 
the loss of the integration constants rj, rj' in the latter cases, since such a linear time dependence 
in if OT ip leads to time dependence in the stress tensor ()2.10p in the non-Ricci flat RN geometries. 

In the present approximation algorithm there are various possible ways of solving the five finite- 
ness conditions ()3.32p . For example, the first condition (I3.32ap is solved by either 

4 = or (4.11a) 
4 = -yAo. (4.11b) 
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FIG. 1: The expectation value (T*^) for a conformally invariant scalar field in the Reissner-Nordstrom 
geometry with Q = 0.99A/. The solid line corresponds to the auxiliary field stress tensor and the dashed 
line to the numerically computed exact one. 

Likewise, independently of this choice, the second condition (|3.32b[l is solved by either 

£i = or (4.12a) 
h = -yAi . (4.12b) 

When the the remaining three conditions are considered, it becomes apparent that either Aq and 
Ai are both zero, or they are both non-zero. The first option, viz. 

4 = Ao = ^1 = Ai = , (minimal) (4.13) 

we term "minimal" because it requires only four integration constants be set to zero in order to 
satisfy all five finiteness conditions (j3.32p . Since this leaves the most integration constants still free 
to adjust, we consider this minimal solution of the finiteness conditions in order to compare the 
auxiliary field algorithm to the direct evaluation of (T^,). The other possibilities for solving the 
conditions (|3.32|) are catalogued in the Appendix. 

The remaining three nontrivial integration constants are ag, ai, and Ci. If the stress tensor is 
finite on the horizon then T\ = T\ there. So once the value of T^* is known, the value of T^q on 
the horizon is fixed by the value of the trace anomaly there. One can compute the components of 
the approximate stress tensor on the horizon from the general expressions ()2.6p . (j2.7p . and (j2.10p 
by using the metric (j3.ip and metric functions (j3.9p along with the expansions (j3.29p and (j3.3ip . 
One finds that the value of r*j on the horizon depends on ai and ci and the value of on the 
horizon depends on oq, ai and ci. 
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FIG. 2: The expectation value (r*^) for a massless spin 1/2 field in the Resisner-Nordstrom geometry with 
Q = 0.99Af. The solid line corresponds to the auxiliary field stress tensor and the dashed line to the 
numerically computed exact one. 
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FIG. 3: The expectation value (T^^) for a conformally invariant scalar field in the Reissner-Nordstrom 
geometry with Q = Q.99M. The solid line corresponds to the auxiliary field stress tensor and the dashed 
line to the numerically computed exact one. 

A numerical code was developed which for a given value of ai chose the values of gq and ci so that 



Tt* and T"" matched the values previously obtained 



from exact numerical computations 



of the stress tensor on the horizon for the spin and spin 1/2 fields. This code then solves the 
equations for the auxiliary fields (p and ^ and computes the analytic approximation ()2.10p for the 
stress tensor for various values of the radial coordinate r. Different values of ai lead to different 
behaviors at large values of r. The goal is to find the value of oi which gives the same energy 
density at large r as the field has if it is in the Hartle-Hawking-Israel state. This is the best that 
one can do to find an approximation for the stress-energy which is finite on the horizon and has 
the correct energy density at large values of r. One could of course take the opposite approach and 
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FIG. 4: The expectation value {T\) for a massless spin 1/2 field in the Resisner-Nordstrom geometry with 
Q — 0.99A/. The solid line corresponds to the auxiliary field stress tensor and the dashed line to the 
numerically computed exact one. 

fix the behavior at large r more accurately, but then the stress-energy would not be regular on the 
event horizon. 

Our results for the spin field are shown in Figures 1 and 3. It can be seen that the approx- 
imation is a reasonably good approximation close to the horizon and at intermediate values of r. 
However, at larger values of r the components of the approximate stress tensor do not approach 
their asymptotic values nearly as quickly as they should. For the spin 1/2 field there were actually 
two values of ai which gave the correct energy density at large r. However, only one is a good 
approximation to the stress-energy tensor near the horizon and that one is shown in Figs. 2 and 4. 

C. ERN Spacetime 

Finally in the ERN case of Q = M the tensor (I2.10p is still able to return a completely finite 
approximation to {TJ') at s = 0, provided that all six conditions in Eqs. (j3.34p and ()3.35p are 
fulfilled. There are three algebraically distinct ways of solving these equations, which are described 
in detail in the Appendix. The simplest possibility which is also minimal in the sense that it 
requires only five integration constants be fixed in order to satisfy all six conditions (j3.34p is the 



31 



3000 
2000 
V 1000 



1 000 



-2000 




12 3 

(r-M)/M 

FIG. 5: The expectation value {T\) for a conformally invariant scalar field in the extreme Reissner-Nordstrom 
geometry. The solid line corresponds to the auxiliary field stress tensor and the dashed line to the numerically 
computed exact one. 
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FIG. 6: The expectation value {T\) for a massless spin 1/2 field in the extreme Reissner-Nordstrom geometry. 
The solid line corresponds to the auxiliary field stress tensor and the dashed line to the numerically computed 
exact one. 
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FIG. 7: The expectation value {T\) for a conformally invariant scalar field in the extreme Reissner- Nordstrom 
geometry. The solid line corresponds to the auxiliary field stress tensor and the dashed line to the numerically 
computed exact one. 
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FIG. 8: The expectation value {T\) for a massless spin 1 /2 field in the extreme Reissner- Nordstrom geometry. 
The solid line corresponds to the auxiliary field stress tensor and the dashed line to the numerically computed 
exact one. 

Despite having two remaining integration constants free, this and the other two solutions of the 
finiteness conditions ()3.34p - ()3.35p in the ERN case are rather restrictive. For example, it turns out 
that for all three of the distinct solutions to (j3.35p . the diagonal components T\,T^j., and T^g on 
the horizon are all fixed (for fixed values of b and b') to the specific finite values, 




rp e 



13M4 M4 
178 , 



13M4 M4 

The correct values of these components on the ERN horizon are {40I 
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FIG. 9: The expectation value {Tg^) for a conformally invariant scalar field in the extreme Reissner- 
Nordstrom geometry. The solid line corresponds to the auxiliary field stress tensor and the dashed line 
to the numerically computed exact one. 
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FIG. 10: The expectation value {Tg^) for a massless spin 1/2 field in the extreme Reissner- Nordstrom 
geometry. The solid line corresponds to the auxiliary field stress tensor and the dashed line to the numerically 
computed exact one. 



The mismatch of the first term in both members of (j4.15p . which is traceless, accounts for the large 
discrepancy of the approximation from the numerical data apparent in Figs. 5-12. This tells us that 
we are certainly lacking some term in Sinv in our minimal approximation based on the anomaly, 
which is needed to give the correct finite coefficient of the stress tensor on the ERN horizon. The 
term ()4.10p is an example of just such a term we have not considered. 
On the other hand, the component depends Hnearly on oq and ci 



1 4484 
^(306856 - 9360ao + 12168ci)6 - ^^^6' 



(4.17) 



585M4' ' 75M4 

which are still free. Hence it is possible to adjust the value of the T"" component on the horizon to 
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FIG. 11: The expectation value (T"") for a conformally invariant scalar field in the extreme Reissner- 
Nordstrom geometry. The solid line corresponds to the auxiliary field stress tensor and the dashed line to 
the numerically computed exact one. 
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FIG. 12: The expectation value (T"^) for a massless spin 1/2 field in the extreme Reissner- Nordstrom 
geometry. The solid line corresponds to the auxiliary field stress tensor and the dashed line to the numerically 
computed exact one. 



any finite value and still have one free parameter left unspecified. This last parameter can be fixed 
by requiring that the behavior of all components of the stress tensor vanish as r — > oo, consistent 
with a zero temperature state in an asymptotically flat spacetime. 

A numerical code was developed which for a given value of ci chose the value of ao so that T"" 
matched the values previously obtained 1^, [3] from exact numerical computations of the stress 
tensor on the horizon for the spin and spin ^ fields. As in the RN case described above, this code 
then solves the equations for the auxiliary fields (j) and ip and computes the analytic approximation 
()2.10p for the stress tensor for various values of the radial coordinate r. 

Our results for the spin and | fields are shown in Figs. 5 -12. Comparison of Eqs. (|4.15p and 
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(j4.16p shows that the T\, T^^, and T^q components are not accurate at the horizon. The plots show 
that this inaccuracy also occurs at intermediate values of r. Since it is possible to choose parameters 
to fit the correct value of T"" on the horizon, the approximation works better for that component. 
Because of the more severe horizon divergences possible in the ERN case, the finiteness conditions 
(|3.34p - ()3.35p eliminating them are restrictive enough to lead to stress tensors which are not very 
good approximations to the numerical results obtained by the direct method. One possibility is to 
consider auxiliary field solutions for which there is a logarithmic divergence in T"" but no power 
law divergence. In this case it is possible to obtain the correct values of Tj*, Tr^, and T^g near the 
horizon. However we have found that the approximation remains poor at intermediate and large 
values of r and of course it also gives a divergence in T"" on the horizon. 

A more promising approach is to add additional terms in Sinv, as for example (j4.10p . in an 
attempt to find a better approximation that is finite on the horizon. Although a modification such 
as this is likely to introduce enough new parameters to allow for significant improvement of the 
comparisons with the numerical results in Figs. ([5] |12p . we do not pursue this possibility here, 
preferring to put the minimal two field anomaly action and stress tensor to its most stringent test. 
We also did not add local terms to the effective action, such as CahcdC"'^'^'^ , which certainly should 
be present in general, and which give contributions to (T'^) of the same order as the anomaly action 
in regular states (where all contributions are of order and small for M ^ Mpi). While the 

approximation based solely on the anomaly action (j2.4p - ()2.5p itself is not very numerically accurate 
for ERN spacetimes, it is nevertheless worth emphasizing that it is the first finite approximation that 
has been obtained for massless quantized fields in these spacetimes, and therefore captures some 
element of the exact effective action lacking in previous approaches, even without the improvements 
possible with the addition of local or Weyl invariant terms. 

V. CONCLUSIONS 

The effective action associated with the trace anomaly provides a general algorithm for approx- 
imating the expectation value of the stress tensor of conformal matter fields in arbitrary curved 
spacetimes. It successfully classifies the leading and sub leading divergent behaviors of the quan- 
tum stress tensor components in the vicinity of all spherically symmetric event horizons, and the 
conformal properties of horizons. These behaviors follow from an analysis of the allowed singu- 
lar behaviors of the auxiliary fields as the event horizon is approached, which are determined by 
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solutions of the conformally invariant differential operator (j2.9p . 

A numerical solution of the auxiliary field equations is necessary to construct the stress tensor 
at points arbitrarily distant from the horizon, except for the uncharged Schwarzschild case which 
admits a completely analytic solution. Due to the possibility of adding homogeneous solutions 
to the auxiliary field equations of motion (12. Sp . the generic divergent behaviors may be canceled, 
and finite approximations to (T"^) obtained on all charged RN event horizons, including also the 
extremal case of Q = M. 

Although it is possible to constrain the solutions of the auxiliary field equations by the require- 
ment that the stress tensor should be regular on the horizon for all Q, we emphasize that the 
fine tuning of integration constants that is necessary to achieve this suggests that a regular stress 
tensor is very much the non-generic case. This is in accord with the well known fact that to have 
a completely regular static stress-energy tensor in a static black hole spacetime it is necessary for 
the field to be in a Hartle-Hawking-Israel state 3] which is a thermal state at the black hole 
temperature. What may be less well known is that this is also true for the zero temperature ERN 
black hole for which the Hartle-Hawking-Israel and Boulware jsol states coincide. Also the 
fact of having a thermal state at a given temperature (including zero) does not completely specify 
the state. For example one can put boundary conditions on the mode functions at some particular 
value of the radial coordinate r as is done when a spherical mirror surrounding the black hole is 
present 
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The generic diverging behavior of the scalar auxiliary fields and their associated stress tensor 
near the event horizon has a geometric origin in the behavior of the Killing field of the static 
geometry becoming null on the horizon, through (IS.lSp . Thus, possible divergences of the stress 
tensor on the horizon are perfectly consistent with the Equivalence Principle, notwithstanding the 
finiteness of local curvature invariants such as R or RahcdR"'^'^'^ there. 

The comparison of the anomaly induced stress tensor, viewed as an approximation to in 
Eq. (|2.10p . with the direct evaluation of the renormalized expectation value for fields of spin and 
^ shows that it yields a fair approximation for states which are regular on the horizon, becoming 
less accurate quantitatively both away from the horizon and in the case of the degenerate horizon of 
the extreme Q = M case. Although the finite terms may be more accurately obtained by a better 
understanding of the Weyl invariant non-local terms in the quantum effective action, such as (j4.10p . 
the fact that the minimal anomalous terms can give all the possible divergent or regular behaviors 
of the stress tensor in the vicinities of all RN black hole event horizons allows us to regard it as 
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a candidate geometrical action for meaningful backreaction calculations. Since the effective action 
and auxiliary field stress tensor of the anomaly can be computed in principle in any spacetime 
without regard to special symmetries, backreaction calculations are possible with Sanom of (12. 4p 
added to the classical Einstein-Hilbert and matter field actions. Moreover, since the resulting 
system of equations can be treated by classical methods, dynamical backreaction calculations in 
time dependent and even non-spherically symmetric spacetimes undergoing gravitational collapse, 
taking into account one-loop vacuum polarization and particle creation effects, would seem to be 
practically feasible by this approach for the first time. 
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APPENDIX A: GENERAL SOLUTIONS OF THE STRESS TENSOR FINITENESS 

CONDITIONS 

In this Appendix we list the possible ways for solving the horizon finiteness conditions in each 
of the three cases considered in the text. 



1. Schwarzschild Spacetime 

The solution of the regularity conditions in Schwarzschild spacetime can be classified into three 
different groups, depending upon the way one solves (|3.22bl) . 

The first group is special, because one can satisfy the three conditions (13.22P , by specifying only 

two of the free constants, namely q = q' = 2. This is the solution considered in [2]. All of the 
finiteness conditions can be satisfied by specifying only five integration constants viz., 

q = 2 (Ala) 

q' = 2 (Alb) 

2hdH = [hqr]' + bq'r] + h' qrj) - 2(6 + b')cH (Ale) 

2hr]' = [bqr]' + bq'r] + b'qrj) - 2(6 + 6')r? (Aid) 

(bqi + bq'ri + b'qr]) 

CH = V or CH = 26T6^^ ^ ' ^ ' 

The other two groups of solutions use six instead of five constants to satisfy the regularity 
conditions: 

q = 2 (A2a) 

CH = (A2b) 

Coo = Y (A2c) 

dH = (A2d) 

, {bqrj' + bq'r] + b'qr]) 2b' + bq' 
^ = 26 2^" 

, = or ^^ ibW+^Wv + b'qy) ^^2f) 
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and 



q = 2--{q'-2) (A3a) 
2b 

CH = --^dn (A3b) 
20 2b / 20 



- 9 f^'-^' 

b' ( {bqi + bq'r] + ^gT?) - t?(46 + 2b' - bq') \ 
- 2b [ 2b + b'-bq' J ^^^'^ 

V = or ^^jbiv'+^bq^ + b'qy) ^^^^^ 

In the last two groups it is assumed that q' ^ 2. Note that with all of the above set of constants 
the regularity conditions for the stress-energy tensor are satisfied for any luminosity L, which is 
given by 

^ = ^ (bqr]' + bq'7] + b'qr]) . (A4) 

This combination of parameters appears in several places in the above solutions, so that if L is set 
equal to zero, as it is in the Hartle-Hawking-Israel state, considerable simplification of the relations 
above results. 



2. General RN Spacetime 

The simplest possibility for solving all five conditions (j3.32p in the < Q < M cases is the 
minimal one (|4.13|) . It requires the fixing of only four integration constants whose values are given 
in Eq. (|4.13|) . The other two distinct possibilities are: 



^0 

h 
Ai 
ai 
Ao 



2b. 



66' Ao 
e(86Ao - b') 

3 2b 

- - TT 

e b'e 

3(l-e)(-66e + (l-3e)feO 
26(3 - 12e + 20e2) 



3)e + 2Ao(l-3e)] 



(A5a) 
(A5b) 

(A5c) 

(A5d) 

(A5e) 



40 



and 

26 

4 = -yAo (A6a) 

^1 = (A6b) 

. 6(266 + 60Ao , . „ , 

= - .(86A0 + 6O ^^'^^ 
3 46 

ai = - + ^Ao(l-3e) (A6d) 

^ 36^(1 -6)[266(l-2e) + 6^(1 -36)] 
° 26[26e(8e2 + 3) + 6'(3 - 12e + 20e2)] ' ^ ' 

These latter two solutions use five constants and are more restrictive than the first minimal one in 
terms of four integration constants. In particular the values of the components T\,T\ and Tq^ are 
all fixed on the horizon in the latter two solutions. There remains the freedom to adjust T"" on 
the horizon, which depends upon the value of oq. 



3. ERN Spacetime 

In addition to the solution ()4.14p to the regularity conditions (|3.34p - (j3.35p given in the text, two 
other distinct solutions of these conditions are possible, viz. 

a-2 = ~yc_2 (A7a) 
26 

a_i = -yc_i (A7b) 

6062(3 + ci)c_i + 266'(45 + 15ci + 81c_i + 13cic_i) + 26'2(45 + 13c_i) 
" 3066'c_i + 6'2(15 + 13c_i) ^ ' 

306c_i + 6^(15 + 13c_i) 

= 8(246 + 136') ^^^^^ 

35 10626' + 13296'26 + 136'^ ± |6'(246 + 136')| V3240062 + 432066' + 6'2 
^"^ ~ 162063 + 3156626' + 11056'26 ' ^ ^' 



The two solutions differ only in the sign of the square root in the expression for c_i. The 
common feature of these solutions is that only five out seven available integration constants are 
used. Nevertheless, as in the previous RN case, the values of the components T\,T\. and T^g are 
fixed on the horizon for all the finite solutions. There remains only the freedom to adjust T"" 
which depends on and ci, both of which are still free. Below we also list the relevant parameters 
corresponding to the solutions found. 
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For the first solution, (j4.14p used in the text the values on the horizon for a single conformal 
scalar field are 

2757 



907r284M4 



and 



-92250 + 3072ao - 3994ci 



(A8) 



(A9) 



907r284M4 

For the second solution with the positive sign of the square root in ()A7p . and for a single 



conformal scalar field, 



The values on the horizon are 



a_2 ^ -0.005526 

a_i ~ 1.135 

ai ~ -93.83 + 6ci 

c_2 ^ -0.000921 

c_i ~ 0.1892 . 

746 



gOvr^S^M^ 



and 



-465000 + 3072ao + 31070c 



907r284M4 

For the second solution with the negative sign of the square root in 



conformal scalar field. 



The values on the horizon are 



a_2 ^ -3.274 

a_i ~ -18.90 

ai ~ —5.306 + 6ci 

c_2 ^ -0.5457 

c_i ~ -3.151 . 

10320 



907r284M4 



(AlOa) 
(AlOb) 
(AlOc) 
(AlOd) 
(AlOe) 

(All) 

(A12) 
7]), and for a single 

(A13a) 
(Al3b) 
(Al3c) 
(A13d) 
(Al3e) 

(A14) 
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and 

25700 + 3072ao - 22980ci 



907r284M4 



(A15) 
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